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The rainbow schlieren deflectometry has been combined with the computed tomography to obtain three-dimen- 


sional density fields of shock containing free jets and we call the method the schlieren CT. Experiments on the 


schlieren CT have been performed at a nozzle pressure ratio of 4.0 by using an axisymmetric convergent nozzle 


with an inner diameter of 10 mm at the exit where the nozzle was operated at an underexpanded condition. Multi- 


directional rainbow schlieren pictures of an underexpanded sonic jet can be acquired by rotating the nozzle about 


its longitudinal axis in equal angular intervals and the three-dimensional density fields are reconstructed by the 


schlieren CT. The validity of the schlieren CT is verified by a comparison with the density fields reconstructed by 


the Abel inversion method. As a result, it is found that excellent quantitative agreement is reached between the 


three-dimensional jet density fields reconstructed from both methods. 
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Introduction 


To understand the prominent characteristics of com- 
plex shock-containing jet flows, schlieren methods [1] 
are widely used as optical tools because of its simple 
optical arrangement with a high degree of resolution and 
ability to easily observe such structures as jet shear layers, 
a Mach disk, barrel shock waves, slip streams, Prandtl- 
Meyer compression and expansion waves in supersonic 
flows. Also, this method for flow visualization does not 
require the introduction of additives into the flowfield 
and is capable of providing useful qualitative information 
on the variations in fluid density, temperature, and static 
pressure. The conventional schlieren techniques have 
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been extensively employed for qualitative flow visualiza- 
tion. 

On the other hand, the simplest method for visualizing 
jet flow fields with varying refractive index quantitative- 
ly is probably the rainbow schlieren deflectometry [2, 3] 
which has been widely developed by Agrawal and co- 
workers. Kolhe and Agrawal [4] measured the flow cha- 
racteristics of the shock-containing jet issued from a 
thick orifice plate with sharp upstream edges by the 
rainbow schlieren deflectometry. Their study was proba- 
bly the first application of the rainbow schlieren deflec- 
tometry for shock containing jets. Yamamoto et al. [5] 
applied the rainbow schlieren deflectometry for correctly 
expanded supersonic jets from an axisymmetric conver- 
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gent-divergent nozzle with a design Mach number of 1.6. 
They showed that the density data from the rainbow 
schlieren for a shock-free jet are in good quantitative 
agreement with those obtained from their numerical si- 
mulation and Pitot tube measurements. 

Agrawal et al. [6] applied the rainbow schlieren def- 
lectometry with tomography for a three-dimensional 
temperature field in a low speed gas jet. Their study de- 
monstrates that the rainbow schlieren tomography can be 
used for quantitative measurements of density in axi- 
symmetric flows. However, there has been no practical 
application of rainbow schlieren deflectometry for a 
quantitative three-dimensional density field measurement 
in a shock-containing jet. Therefore, the aim of the 
present study is to develop a three-dimensional rainbow 
schlieren system for shock-containing jets. For the pur- 
pose, the conventional rainbow schlieren deflectometry is 
combined with the computed tomography and we call the 
method the schlieren CT in the present paper. The density 
field of the underexpanded sonic jet from an axisymme- 
tric convergent nozzle has been measured by the schlie- 
ren CT. The density field for the same jet flow is also 
displayed from the so-called Abel inversion method 
which is extensively used for the reconstruction of the 
radial density distribution of a cylindrically symmetric 
object. Since this method is very effective for quantita- 
tive density measurements in axisymmetric shock con- 
taining jets [4,5], the density field reconstructed from the 
Abel inversion method is used to confirm that recon- 
structed from the schlieren CT. 


Experimental appratus 


Experiments have been performed in a blow-down 
supersonic wind tunnel with the jet issued in the quies- 
cent laboratory air. A schematic diagram of the experi- 
mental apparatus with a schlieren optical system is 
shown in Fig. 1. The air supplied by a compressor that 
pressurizes the ambient air up to 1 MPa is filtered, dried 
and stored in a high-pressure reservoir consisting of two 
tanks with a total capacity of 2 m°. The high-pressure dry 
air from the reservoir is stagnated in a plenum chamber 
and then discharged into the atmosphere through a test 
nozzle. The total temperature in the plenum chamber is 
equal to the room temperature, and the plenum pressure 
is controlled and maintained constant during the testing 
by a solenoid valve. The test nozzle used in the present 
experiment is shown in Fig. 2. It has an axisymmetric 
wall contour with a 10 mm inner diameter at the exit 
plane and the nozzle wall contour over a range of parts A 
to B is designed by a sinusoidal curve to provide uniform 
parallel flow with respect to the nozzle axis at each loca- 
tion of the A and B. The jet issued from the nozzle is vi- 
sualized by the rainbow schlieren deflectometry with a 
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field of view of 100 mm diameter for a nozzle pressure 
ratio of 4.0, i.e., jet flow in underexpanded condition for 
the present nozzle. 

The rainbow schlieren system consists of rail-mounted 
optical components including a 50 um diameter pinhole, 
two 100 mm diameter, 500 mm focal length achromatic 
lenses, a computer generated 35 mm wide slide with col- 
or gradation in a 1.4 mm wide strip, and a digital camera 
with variable focal length lens. A continuous 250 W met- 
al halide light source connected to a 50 um diameter fiber 
optic cable provides the light input at the pinhole through 
a 16.56 mm focal length objective lens. The camera out- 
put in the RGB format is digitized by a personal comput- 
er with 24 bit color frame grabber. The rainbow filter 
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Fig.2 Schematic drawing of test nozzle 
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used in the present experiments is shown in Fig. 3. The 
filter was fabricated in computer software and then 
printed digitally on a high resolution 23 mm color film 
recorder. It has continuous hue variation from hue = 0 to 
314 deg in a 1.4 mm wide strip and the origin y = 0 cor- 
responds to Hue = 180 deg. The characteristics of the 
rainbow filter are performed by traversing the filter au- 
tomatically in intervals of 0.1 mm in the y direction at the 
schlieren cut-off plane before starting experiments. The 
calibration result is shown as open symbols in Fig. 4. The 
abscissa is Hue and the ordinate is the transverse ray dis- 
placement d from the x axis at the cut-off plane. The solid 
line indicates a least-squares regression line of the expe- 
rimental data using fifth degree polynomials. 

As shown in Fig. 5, the test nozzle is installed inside a 
pulley with a gear ratio of 1:4 and it can be rotated about 
the center axis of the nozzle by a stepping motor with a 
gear ratio of 1:36 connected through a timing belt. In the 
present experiment, rainbow schlieren pictures are ac- 
quired over a range of nozzle angular angles from 0 = 0 
deg to -180 deg by rotating the nozzle about its longitu- 
dinal axis (z axis) in equal angular intervals of -10 deg. 
These 19 schlieren pictures are used for reconstruction of 
the jet density field by the Abel inversion and schlieren 
CT methods. The jet density profiles obtained from both 
methods are compared with each other. 
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Fig. 4 Calibration curve of rainbow filter 
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Rainbow schlieren visualization 


A typical rainbow schlieren picture of underexpanded 
sonic jet issued from an axisymmetric convergent nozzle 
is shown in Fig. 6 with the flow from left to right. The 
ratio Pos/P» of the plenum pressure pos to back pressure 
P» is held constant at 4.0. As shown on the upper right 
corner in the picture of Fig. 6, the rainbow filter is placed 
horizontal with respect to the z axis at the schlieren 
cut-off plane and the background hue corresponds to the 
color of the location indicated as the right pointing arrow 
for the filter image. The schlieren picture has a high res- 
olution of 55 pixels/mm and it corresponds to a spatial 
resolution of around 20 um for the present schlieren set 
up. The schlieren picture of Fig. 6 shows the archetypal 
structure of underexpanded jet including the expansion 
fans, barrel shock, Mach disk, reflected shock, and slip 
stream. Also, a significant color gradation along the radi- 
al direction can be clearly seen near the jet boundaries of 
y =+5 mm and - 5 mm. Background hues can be seen at 
the jet centerline and in the regions radially away from 
the jet boundaries where the transverse component of the 
density gradient is intrinsically zero. Also, the hues cor- 
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Fig. 6 Typical rainbow schlieren picture of underexpanded 
sonic jet for po;/pp=4.0 
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responding to the upper and lower jet boundaries are dif- 
ferent from each other and they exist in the lower and 
upper regions from the background hue in the filter im- 
age shown in Fig. 6, respectively, because an axisymme- 
tric rainbow filter shown in Fig. 3 is used in the present 
schlieren system. Therefore, the light rays passing 
through the upper and lower jet boundaries have down- 
ward and upward angular deflections toward the opposite 
boundaries, respectively. 


Density contour plots 


The density field corresponding to the schlieren pic- 
ture of Fig. 6 can be estimated from two reconstruction 
processes. One is the analytical method based upon the 
Abel inversion and it is extremely effective for the re- 
construction of the radial density distribution of a cylin- 
drically symmetric jet. The other is that based upon the 
schlieren CT and it is valid even for that of an axisym- 
metric jet. The detail descriptions for the reconstruction 
of density fields based upon the Abel inversion method 
and the rainbow schlieren combined with the computed 
tomography are reported in the paper of Kolhe and 
Agrawal [4] and that of Agrawal et al. [6], respectively. 

The density contour plots reconstructed from both 
methods are compared with each other as shown in Figs. 
7(a) and 7(b) where the density contour plot of Fig. 7(a) 
is computed as the mean of the three density fields esti- 
mated from the Abel inversion using the horizontal 
schlieren images recorded at nozzle rotational angles of 0 
= 0, -60 and -120 deg, because, the reconstruction pro- 
cess based upon the Abel inversion method is relatively 
sensitive to noise and errors in the determination of the 
jet centerline [3]. Also, Fig. 7 is presented with the color 
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Fig. 7 Axial density contour plots for po,/pp = 4.0 


representation of density field and the solid lines with a 
numeral in the contour show isopycnic. 

These contour plots indicate good qualitative and 
quantitative indications of the shape of the various fea- 
tures of the shock containing jet structure. The density 
contour of Fig. 7(b) exhibits nearly axial symmetry about 
the jet centerline. From a comparison between Figs. 7(a) 
and 7(b), it is found that the global shapes of each shock 
in the shock cell structure are almost identical to each 
other. In addition, the axial distances of the Mach disk 
from the nozzle exit plane are nearly the same for both 
density fields. 


Density distributions along axial and radial directions 


A comparison of density profiles along the jet center- 
line obtained by the Abel inversion and schlieren CT for 
an underexpanded sonic jet is given in Fig. 8 where den- 
sity data near the nozzle exit are excluded from Fig. 8, 
because in areas very close to the nozzle exit, reflections 
of the schlieren light rays by the exit wall surface lead to 
larger errors in the density field or yield regions where 
density data cannot be precisely obtained. The abscissa of 
Fig. 8 is the downstream distance z from the nozzle exit 
plane and the ordinate shows the density. Also, the blue 
and red lines in Fig. 8 indicate the densities reconstructed 
by the Abel inversion and schlieren CT methods, respec- 
tively, and the dashed line parallel to the abscissa is that 
in ambient air. 

For a nozzle pressure ratio of p,/py = 4.0 in the 
present experiment, the density p, at the nozzle exit plane 
calculated based upon the assumption of the one-dimen- 
sional isentropic flow from the nozzle inlet to the exit 
becomes a value of 3.04 kg/m’ and it is shown as the 
left-pointing arrow on the vertical axis in Fig. 8. There- 
fore, the densities inside the nozzle are always higher 
than the density p, = 1.20 kg/m’ in surrounding air. Since 
the static pressure at the nozzle exit plane is greater than 
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Fig.8 Axial density profiles for x = y = 0 mm 
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the back pressure, the freestream Mach numbers at the 
nozzle exit plane are almost unity and Prandtl-Meyer 
expansion fans occur from the lip of the nozzle exit. A 
conventional shock cell structure appearing repeatedly in 
the jet flowfield can be seen in Fig. 8. Also, the shock 
amplitudes in the shock-cell structure decrease gradually 
toward the downstream direction and there is very fa- 
vorable overall agreement between both density profiles 
over the full length of the region measured. Both the 
shock-cell intervals and density amplitudes agree ex- 
tremely well. 

In order to observe the three dimensional structure of a 
shock-containing jet, the cross-sectional density contour 
plot at any axial location in the jet flowfield was obtained 
from the schlieren CT using the horizontal rainbow 
schlieren images taken over 19 viewing angles from 0 = 
0 to -180 deg by rotating the nozzle about the z axis in 
equal angular intervals of -10 deg. The result for the jet 
density field at an axial location of z = 10 mm is depicted 
in Fig. 9 with the color bar showing density ranges. This 
axial location exists in the region just prior to the Mach 
disk. 

The density contour plot of Fig. 9 exhibits nearly 
symmetry about the center of the contour. However, the 
axis of symmetry does not coincide with x = y = 0. It 
appears that the direction of the optical axis of the rain- 
bow schlieren system in the present experiment may have 
deviated slightly from that of a vector along x axis of Fig. 
1. Also, the isopycnic line with p = 1.6 kg/m’ around 4 
mm away from the center of the contour shows a ring 
with some indentations. It is deemed responsible for the 
presence of shocklets in the flow observed in the study of 
Fourguette et al. [7]. However, Krothapalli et al. showed 
that it results from some imperfections of the nozzle in- 
ner surface [8]. Such three-dimensional flow features 
cannot be easily obtained from the line-of-sight approach 
such as conventional schlieren and shadowgraph tech- 
niques. 
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Fig.9 Cross-sectional density contour plot at z = 10 mm 
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The radial density profiles reconstructed by the Abel 
inversion and schlieren CT are shown in Fig. 10 as the 
blue and red lines, respectively. The profile corresponds 
to the densities obtained along y axis at x = 0 in Fig. 9. As 
can be seen in Fig. 10, the jet has lower densities com- 
pared with the ambient one over a range of radial dis- 
tances from r = 0 mm to around 2 mm, which is due to 
the isentropic expansion from the plenum chamber up- 
stream of the nozzle to just upstream of the Mach disk. 
Also, from the density profile by the schlieren CT, the 
density field just upstream of the Mach disk is almost 
axisymmetric with respect to the horizontal axis includ- 
ing y = 0. Excellent quantitative agreement is reached 
between the density profiles reconstructed from both 
analytical methods. 
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Fig. 10 Radial density profiles at x = 0 mm and z = 10 mm 
Concluding Remarks 


Three-dimensional density fields of the shock con- 
taining free jet issued from an axisymmetric convergent 
nozzle with an inner diameter of 10 mm at the exit were 
quantitatively obtained by the rainbow schlieren deflec- 
tometry combined with the computed tomography and 
we call the method the schlieren CT. As a result, a densi- 
ty contour plot at any jet cross-section can be efficiently 
inferred from the schlieren CT and a density field in- 
cluding the jet axis quantitatively shows the representa- 
tive shock-cell structure of repeated compression and 
expansion waves by continuous color gradation. Axial 
and radial density profiles by the schlieren CT were 
compared with those by the Abel inversion method. The 
density profiles obtained from both methods are quantita- 
tively identical to each other. 
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